clear all
close all
clc


% D0=xlsread('ABAQUS.xlsx');
% D1=load('FssiCAS/TimeHistory');
D1=load('TimeHistory');

St=12;
FS1=16;
subplot(1,2,1);
hold on
plot(-(D1(St:end,10)+D1(St:end,11)+D1(St:end,12))/3000,(D1(St:end,10)-D1(St:end,12))/1000,'color','m','linewidth',2,'LineStyle','--');
% plot(D0(:,3),D0(:,4),'color','r');
% set(gca,'xlim',[0 240]);
% set(gca,'xTick',0:40:240);
% set(gca,'xticklabel',[]);
% set(gca,'ylim',[0 200]);
% set(gca,'yTick',0:40:240);
% set(gca,'yticklabel',[]);
set(gca,'box','off','Fontname','Times New Roman','Fontsize',FS1);
xlabel('p (kPa)','Fontname','Times New Roman','Fontsize',16);
ylabel('q (kPa)','Fontname','Times New Roman','Fontsize',16);


subplot(1,2,2);
hold on
plot(-D1(St:end,6),(D1(St:end,10)-D1(St:end,12))/1000,'color','m','linewidth',2,'LineStyle','--');
% plot(-D0(:,1),D0(:,4),'color','r');
% set(gca,'xlim',[0 0.15]);
% set(gca,'xTick',0:0.03:0.15);
% set(gca,'xticklabel',[]);
% set(gca,'ylim',[0 150]);
% set(gca,'yTick',0:30:240);
% set(gca,'yticklabel',[]);
set(gca,'box','off','Fontname','Times New Roman','Fontsize',FS1);
set(gcf,'Unit','Centimeters','position',[10 10 25 15],'color','w');
xlabel('\epsilon','Fontname','Times New Roman','Fontsize',16);
ylabel('q (kPa)','Fontname','Times New Roman','Fontsize',16);



